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ABSTRACT 

Alfven waves may be difBcult to excite at the photosphere due to low ionization fraction and 
suffer near-total reflection at the transition region (TR). Yet they are ubiquitous in the corona and 
heliosphere. To overcome these difficulties, we show that they may instead be generated high in 
the chromosphere by conversion from reflecting fast magnetohydrodynamic waves, and that Alfvenic 
transition region reflection is greatly reduced if the fast reflection point is within a few scale heights 
of the TR. The influence of mode conversion on the phase of the reflected fast wave is also explored. 
This phase can potentially be misinterpreted as a travel speed perturbation, with implications for the 
practical seismic probing of active regions. 
Subject headings: Sun, oscillations; sunspots; Magnetohydrodynamics (MHD) 
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1. MOTIVATION 



In recent years, Alfven waves have been shown to be ubiquitous in the solar corona fTomczyk et al.|[2007 Mcintosh 



et al. 2011 1 though these may more properly be classed as magnetohydrodynamic (MH D) kink waves wh ere coronal 



magnetic structuring supports such tube waves 
case, they are essentially transverse and incompressive 



20081. In either 



( [Edwin fc Roberts||1983[ [Van Doorsselaere et al 

issive, and may b e reterrea to as "Alfvenic" . Alfven waves have also 

been detected in situ in the solar wind ( Belcher fc Davis Jr||1971 1, of which they are postulated to be crucial drivers 
( [Hollweg 2006). Furthermore, Alfven waves are inferred in t heupper chromosphere based on Hinode Solar Optical 
itelescope (SOI') observed transverse oscillations of spicules dDe Pontieu et al.||200 7j). Althou gh long-period Alfven 

It 



Hollweg"2006 ) , coronal observations 
-4 mHz apparently associated with 



waves (hours) dominate observations in the solar wind ( BelcheTfc Davis Jr| 

reveal oscillations largely in the 100-500 s range, with a discernible peak at 

p-modes. In this paper we shall be concerned only with these shorter period waves, as our main purpose is to explore 

the extent to which coronal Alfven waves may be extensions of the Sun's p-mode wave field. The source of long-period 

waves must be sought elsewhere. 

Traditionally, it has been thought that direct generation at the solar photosphere by granular buffeting is responsible 
for the Alfven waves observed further out in the solar atmosphere and beyond (e.g. Cranmer fc van Ballcgooijcn 2 005 ■). 
Howe ver, this does not take account of the i nefficiency of photos pheric MHD Alfv6n generation (p arkcr 1991; C ollins] 
19921, nor of low ionization fraction effects (Vranjes et al. 2008 though see Tsap, Stepanov, fc "Ko 
contrary view)Fl 

Another impediment to Al fven waves reaching the corona is the strong reflection they suffer at the chromosphere- 
corona tran sitio n regio n (TR; Uchida fc Sakurai 1975 1. This effect can even produce a trapped- wave resonance structure 
(jSchwartz et al.||1984|), although as pointed out by Cranmer fc van Ballegooijen| ( |2005[) it is somewhat detuned by 

so completely absent in the case ol: an open lower boundary as 



pylova, |2011| for a 



although as pointed out by 

non-Isothermal stratification. The resonances are a' 



adopted here. Cranmer fc van Ballegooijen quote typical Alfven reflection coefficients of around 95%, and we concur in 
general though perhaps we w ould put that figure even higher at around 98% in the cases discussed here. Nevertheless, 
Cranmer & van Ballegooijen conclude that this may still be enough to sup ply chromospheric and coronal energy losses 
and to power the solar wind, contrary to the view of Rosner et al. ( 1986 K_ 



Photospheric Doppler and Zeeman observations at two heights ( Ulrich 1996 ) suggest that indeed there may be 
sufficient Alfven flux in the 5-minute band to supply coronal losses in the quiet Sun. Indirect arguments suggest that 
even the order-of-magnitude greater energy requirements of active regions may be consistent with these observations, 
though this conclusion relies on extrapolation. The velocity data are strongly indicative of outgoing Alfven waves, 
though the noisier magnetic data may suggest the addition of a reflected (downgoing) component, which might be 
expected given the high reflectivity of the transition region. 

Irrespective of this, the TR is undoubtedly a considerable hurdle for Alfven waves, and it remains to be seen whether 
the observed coronal Alfven amplitudes are consistent with waves originating in the photosphere alone in both quiet 
and active regions. Is there some way that the Alfven transmission coefficient can be increased'ijj 

We postulate and explore a novel possibility: that the Alfven waves have not propagated as transverse waves (Alfven 
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'- Of course, Alfven waves may be produced in isolated structures in the low atmosphere. For example, [Jess et al.| l |20C)9| l identify torsional 
Alfven waves in a bright-point group of kilogauss strength, though the height of excitation is not clear from this study, and ionization 
fraction may be enhanced compared with quiet Sun. ^^_^^^^_^_ 

^ The situation in coronal loops may be different, with resonances allowing greatly enhanced transmission | |Hollweg|l99l| . We focus on 
open field only though. 
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or kink) from the photosphere, but that instead they have been generated by local mode conyersion from reflecting 
fast waves high in the chromosph ere ( Melrose|1977 Melrose fc Simpson|1977 Cally fc Goossens|2008 Cally 
] Khomenko fc Cally| 2012), where low ionization fraction is not an issue tor the 3-6 mHz trequenc; 
It wif 



2011 

most Interest here 



& Hansen 

ency range of 

be shown that this greatly increases Alfven penetration of the TR, thereby opening up a 



potentially more fruitful source of coronal Alfven waves. 

The progenitor fast waves bear the signature of this process not only in their partial loss of energy to the Alfven 
wave s but also in their phase as they return to the interior. Local heliosei smology uses phase to infer wave travel times 
(e.g., 
effect 



Duvall Jr et al.|["l993 Braun fc Lindsey||2000 Gizon fc Birch|[2005 1 , so it is of interest to gauge this anomalous 
I'his is discussed in Section 14^ 



2. FAST-TO- ALFVEN CONVERSION: REPRISE OF PAPER I 

The solar interior is populated by p-modes, global modes of oscillation excited by convection near the surface and 
trapped in a frequency-dependent cavity between an acoustic cutoff at the top and the Lamb depth at the bottom. 
At frequencies above about 5 mHz the acoustic cutoff barrier is breached and the waves can escape into the solar 
atmosphere. Strong inclined magnetic fields, as found in sunspots and other smaller mag netic field concentrations, can 
reduce the effective cutoff to allow even lower frequency acoustic waves to escape ( [Jefferies et al.|[2006j . 

Now, in magnetic regions the Alfven speed a increases rapidly with height due to density stratification. So at 
some point, inevitably, the Alfven speed surpasses the sound speed c. As the waves pass through the Alfve n-acoustic 
equipartition layer where a and c coincide, mode conversion can occur ( Schunker fc Cally]|2006 Cally|[2007 ), splitting 
the waves into fast (i.e., predominantly magnetic) and slow parts. As the slow wave propagates higher it becomes 
progressively more field-aligned and more acoustic in nature. Conventionally, the resultant slow waves are said to be 
caused by transmission and the fast modes by conversion. The amount of transmission/conversion here is dependent 
on the attack angle a, the angle between the wavevector and magnetic field direction. A small attack angle allows 
greater transmission to the slow wave and correspondingly a larger attack angle produces greater conversion to the 
fast wave. In an exact isothermal model conversion to the fast wave can be near total at large a (Hansen fc Cally 



2009f. 

Here we assume that this process has taken place, and follow the fast wave as it progresses to ever higher Alfven 
speeds a ^ c. In this regime we adopt the simplifications afforded by the cold plasma approximation and set c = 0, 
or equivalently /3 = where t he plasma-^g represent s the ratio of gas to magnetic pressures. 

In Paper I of this sequence ( Cally fc Hansen|2011 ) we investigated fast-to- Alfven mode conversion in a simple /3 = 
plasma with uniform inclined magnetic field Bo{cos9,0,sm9) and exponentially decreasing density p ex e^^''* in the 
direction x of inhomogeneity, where h is the scale height. This may be thought of as the vertical direction in a plane 
stratified model of the solar atmosphere in which density decreases exponentially with height. A fast magnetohydro- 
dynamic wave is injected from x — > — oo, propagates to the right (positive x direction) and reflects (classically) at 
w^ = a^(fcy -|- kl), where a{x) ex e^^'^'^ is the Alfven speed. An exp[i{kyy + kzZ — ojt)] time and transverse space 
dependence is assumed. The orientation of the wavevector in y-z space is arbitrary. 

If ky =/= 0, the fast wave partially converts to Alfven waves, typically around and beyond the fast wave reflection 



point. Introducing the dimensionless transverse wavenumber k — \/k'1 + k^, where Ky = kyh ~ Ksini 

K cos (f), it was found that very significant fast-to- Alfven conversion occurs in various regions of K-6-(f> 
Specifically, it was found that: 



b and k^ — k^h — 
parameter space. 



1. The thickness of the conversion region depends sensitively on k. For k ^ 1 it is several scale heights thick, for 
K ^ 0.2 it extends some 20 scale heights beyond the reflection point, but for k ~ 5 it is less than one scale height 
wide. For most solar atmospheric waves of interest in the chromosphere, we might expect k to be considerably 
less than 1, indicating that the Alfven conversion region can effectively flU the chromosphere. 



For (f> < 90°, conversion is to outgoing (upward) Alfven waves, but for cj) > 90° conversion predominantly 
produces backward (downward) propagating Alfven waves. This was to be expected based on a consideration of 
whether the fa st wave aligns with the magnetic field before or after its reflection, consistent with the flnding of 
Melrose ( 1977 ) that fast-to- Alfven coupling is strongest at small attack angle. 



These predic tions were amply conflrmed by numerical simulations in a more realistic model sunspot by [Khomenko fc 
Caily] (|2012l) 



This simple one-layer model crudely represents the solar chromosphere, and takes no account of the transition region 
(TR) or corona above. Nevertheless, it allowed us to better understand the local conversion mechanism that goes on 
there without complications arising elsewhere. In this paper though, we extend the model to consist of two exponential 
layers, representing the chromosphere a^ ex e^^^ and the corona a^ oc e^^'^""' , with /icor ^ h and a discontinuous jump 
between them representing the TR. This will allow us to explore the consequences of Alfven reflection off the TR and 
truncation of the extended conversion region by the TR. 

3. MODEL, EQUATIONS, AND NUMERICAL APPROACH 

In sunspot umbrae the Alfven-acoustic equipartition level a = c i s typically situated a few hundred kilometre s 



below the photosphere and at about the photosphere in the penumbra (Mathew et al.|2(304 Borrero fc Ichimoto|2011 ) 
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Around active regions and in mor e diffuse magnetic environments equ ipartition is normally identified with the magnetic 
canopy of the low chromosphere ( Bogdan et al.|[2003 Finsterle et al.j 2004 )pl As explained in Section^ seismic waves 
from the solar interior split at a w c into slow (acoustic) waves and fast (magnetic) waves that may propagate into 
the atmosphere above. We are not concerned with the slow wave here, but the fast wave's fate is very interesting as 
it can suffer a further mode conversion to the Alfven wave (Paper I). 

Seismic waves below about 3~4 niHz (depending on field strength) may not reach the a — c equipartition level before 
reflecting. A precondition of the study here is that they have passed through a — c, and therefore partially mode 
converted to magnetically dominated fast waves in a ^ c. 

With this in mind, we now move well above the a ^ c region and adopt a uniform-magnetic-ficld two-isothermal-layer 
Alfven speed profile 



A^ {h^or-/h)e<-'=-''^^^''''°' X > XT, 



(1) 



crudely representing the low-/3 chromosphere (x < xt) and the corona {x > xt) separated by a transition region 



discontinuity. Here x is the vertical coordinate, increasing with height in the solar atmosphere, h 



'ha and 



^cor = c^o-c/ld ^'''6 the density scale heights in the two regions, c and Ccor are the respective sound speeds, 7 is the 
adiabatic index, and g is the gravitational acceleration. The factor A is the Alfven speed at the top of the chromosphere. 
The exponential increase in a{x) results solely from an exponentially decreasing density p{x). The density jump across 
the TR is of course equal to the inverse of the temperature jump, which explains the hcor/h factor in the second line 
of Equation M. Typically, /icor 3> h. 

The crude uniform field model is more applicable to large-scale strong magnetic field regions such as sunspots than to 
network or other structures where discrete flux tubes at the photosphere quickly expand with height in the atmosphere 
until they abut and thereafter become more uniform. Especially in this latter case, the exponential Alfven speed 
increase supposed here will be partially ameliorated by the geometric field spread. Nevertheless, the model represents 
a useful testbed on which to discern general principles. 

To focus on the fast-to- Alfven conversion process in the upper solar atmosphere it is convenient to ignore acoustic 
waves altogether by neglecting the sound speed c compared to the Alfven speed a in the wave equations, the so-called 
cold plasma or zero-/3 approximation. The underlying equilibrium stra tification is unaffected. This d evice is used 
commonly in the study of MHD waves in the upper solar atmosphere (Nakariakov & Verwichte 2005). There are 
innumerable papers in the literature that have studied MHD waves in the cold plasma regime, representing the Sun's 
outer atmosphere, and have adopted nonu niform density distri butio ns to create Alf ven speed inhomogeneities, e.g., 
Davila ( 1987|) , Ruderman & Roberts ( 2002| ) ,|Pascoe et al. (2 010|), a nd Hollweg (1990) amongst many others. This last 
reterence is closest to our study, and tnat 01 (Jally & Andries 12010), in that it discusses a fast MHD wave injected into 
a region of increasing Alfven speed in which it reflects but also transfers energy to an Alfven resonance. The density 
profile there is linear though instead of our expon ential. The close corr e spond ence of the results of Paper I, which uses 
the /3 = approximation, and the simulations of Khomenko & Cally (2012), which do not, supports the utility and 
validity of the cold plasma model in this instance. 



We also assume we are not in the low- frequency regime of atmospheric gravity waves ( Straus et al.|[2008 ). Gravity- 
wave-to- Alfven conversion for very low freque ncy waves (^ 1 mHz, co mparable to the Brunt- V aisala frequency N) has 
been extensively treated elsewhere ( _Newington fc Cally 2010 2011), and is found to be viable only in very inclined 
magnetic fields. We therefore assume w^ ^ N'"' and that buoyancy plays an insignificant role. 

With all this in place then, we adopt the linearized equation governing the oscillations of a cold plasma (c.f.. Equation 
(1) of Paper I) 



9n 



-VpX, 



(2) 



where ^{x) 
subscripts 



^x^x + £.v^v + ^zf^z ~ S,±e± + ^yGy is the plasma displacement and X = ^'^ is the dilatation. The 
and '-l' denote the parallel and perpendicular (in the x-z plane) directions to the uniform magnetic field 

-By Oy. 



Bo = Bo{cos9, 0, sin0) respectively. The 'p' refers to the full perpendicular component to Bq, i.e., Vp = e_L i9_l- 
Equation ^ neatly expresses the role of the fast wave (represented by x) as a source of Alfven waves (the term in 
brackets on the left hand side is the pure Alfven operator). 

Frobenius and WKB solutions to this equation with exponential Alfven speed were developed in Paper I for a; — >■ 00 
and X -^ —00 respectively, and they are again used here. Indeed, the Frobenius solution is utilized throughout the 
corona. The imposed boundary conditions are that there are no incoming Alfven waves from x = ±c», and that the 
exponentially decaying fast wave solution is selected at x ~ +00 with the exponentially increasing solution deprecated. 
Numerical integration from the WKB region to a suitable matching point in the chromosphere and from the TR also 
to this matching point completes the solution. Continuity of $ and d^/dx is applie d across xt- 

It is worth reiterating the results from the Appendices of Cally & Goossens ( 2008 ) and Paper I that the point a; — >■ c» 



is a regular singular point of the wave equations for the coupled fast and Alfven waves. In the case ky = where 
the two wave types decouple, exact solutions in closed form are available. For the Alfven wave these are in terms of 



r(i) 



r(2) 



the Hankel functions Hq and Hq , representing respectively downgoing and upgoing waves (see Equation (|3|) and 

^ Usually such discussions are couched in terms of the plasma beta equals unity layer, where /3 = Pgas/Pmag is the ratio of the gas to 
magnetic pressure. However, = 1c? j^a^ , where 7 is the ratio of specific heats, so the distinction between /3 = 1 and a = c is almost 
immaterial. 
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Appendix [A| for full details). In the general case addressed here {ky ^ 0), no such closed form solutions are available, 
but full Frobenius series solutions were constructed in Appendix A. 2 of Paper I. These allow us to mat c h our numerical 
solutions to the pure outgoing Alfven wave at a; — >■ +00. We note that Cranmer & van Ballegooijen (2005) find that 
reflection of short period Alfven waves is weak in the corona proper, supporting our selection of a radiation boundary 
condition. 

Assuming the incoming fast wave from a; = — 00 carries unit wave-energy flux jF^ = 1, we calculate the outgoing 
fast flux jF^ at a; = —00, the outgoing Alfven flux J^+ at a; = +00, and the outgoing Alfven flux J^^ at x = — cx). Of 
course, Si'^ + J^^ + jT^ — J-^ — 1, which provides a non-trivial test on our numerics. 

Determining the relative proportions of the three outgoing modes is our prime goal. In the single layer model of 
Paper I, any Alfven wave launched in the positive x direction from the conversion region continued outward to yield 
^+ , whereas an initially backward-propagating Alfven wave excited by the reflected fast wave gave A" . Here though 
the picture is complicated by Alfven reflection from the TR, leading us to expect a greatly diminished J1+. Whether 
S\~ increases or decreases is not clear a priori: TR reflected j^+ should enhance Jl~, but truncation of the well-spread 
chromospheric conversion region will diminish it. Numerical solution is required to assess which of these two competing 
processes prevails. 

The full parameter set describing our system is dimensionless wavenumber k, magnetic field inclination 9, wavevector 
polarization 0, transition region position xt, and coronal scale height /icor- The chromospheric scale height h is 
arbitrarily scaled to unity. The frequency w and Alfven speed scaling A can be absorbed into the position of fast wave 
reflection that occurs at Loh/a — k. The zero of x is fixed so that uj'^h^ /a^ — exp[— x/ft.] in the chromosphere, so fast 
wave refiection occurs at xr — —hlnK^ (provided this is less than xx)- We investigate how the various energy fiuxes 
depend on each of the free parameters. 

4. RESULTS 

4.1. Transmission of the Pure Alfven Wave 

Before embarking on a study of fast-to-Alfven mode conversion, it is useful to first calculate the transmission 
coefficient ^ "T ^ 1 of a pure Alfven wave incident on the transition region from the chromosphere. Here 'T 
represents the fraction of wave energy fiux in the incident Alfven wave that passes into and through the corona. The 
remainder is refiected at the TR. The pure wave, with ^y polarization only {(jj — 0), is completely decoupled from and 
therefore not excited by the fast wave, but it does give us a good measure of how strongly the TR reflects Alfven 
waves. The exact solutions of the Alfven wave equation d^^y/dt^ — c? d'^^yjds^ (where s is distance along a field line) 
in an atmosphere with exponential Alfven speed are available in terms of Hankel functions: 



iy = ^^,^e-'^^A--T)t.nejjii,2) f^seci 



(3) 



with H^ representing a wave propagating upward and H^ corresponding to downward travel. Si 2 are arbitrary 
amplitudes with dimensions of length. As explained in Appendix [X] the Hankel functions represent unidirectional 
waves. This is consistent with the total Poynting flux in the x-direction: 



-F. - 



ujBq cos ( 
n h ^ 



F2| 



Fll 



(4) 



r(2) 



Matching to an Hq solution (i.e., outgoing only) with the coronal scale height and letting ft, = 1 in the chromosphere, 
the transmission coefficient is 



T=l- 



r(2) 



1/2 



W^' 2K hcLr sec H\' {2K sec 



r(2) 



-1/2 „(2) 



(2) 



H^q' {2K sec e)H\> 2K KLr sec 6 



1/2. 



hcoV^ hI)^\2K aec0)H'l^^ hKhlihecO) - h'i^\2K sec0)H^^^ (2Khlihec0 



(5) 



where K = ujh/a[xTp) = cxp[— a:T/2/i] is the dimensionless Alfvenic field-aligned wavenumber at the top of the 
chromosphere. Note that, for fixed /icor, '2' depends only on K sec 0. T increases monotonically from zero at K sec — 
towards 

as Ksec9 — >■ 00 (see Fig. II]). For /icor = 100, "I^nax = 0.3306. Only for /icor — 1 (no transition region) does l^iax reach 



For comparison, the pure Alfven transmission coefficient for the case of a uniform rather than exponential corona is 



T = 1 



r(2 



Wq> {2K sec 6) + i hJotHY-' {2K sec 



l/2rr(2)/ 



H''^\2Ksec( 



ihU?H[^\2Ksect 



(7) 



where now /icor represents only the jump in a^ across the TR and not the coronal scale height (which is infinite). 
It results directly from (|5| in the limit 2KhcoT sec^ ^ 1 where the relevant Hankel functions may be replaced by 
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Figure 1. Left: Pure Alfven transmission coefficient CT (see Equation (Isl) against Kscc9 for hcor = 20, 50, 100, 150, and 200 (full curves, 
top to bottom). In the limit hcor —^ 1, the TR vanishes and transmission becomes total. The dashed curves represent the transmission 
coefficient for a uniform corona, as given by Equation (ItI. Right: Contour plot of transmission coefficient T against transition region 
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position XT and magnetic field inclination 9 for the case 



100. The contour heights are as labelled. 
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Figure 2. The forwrard Alfven conversion coefficient ^+ (full curve), reverse Alfven conversion coefficient Si (dashed curve) and reverse 
fast conversion coefficient ^~ (dotted curve), as a function of hcor for k = 0.2 (i.e., xji = — lnO.2^ = 3.219), 9 = 30°, xx = S, (f> = 70° (left 
panel) and </> = 110° (right panel). Both variables are plotted on logarithmic scales. The heavy long-dashed curve represents T as given 
by Eq. kh. The chained lines represent S^T, where J^(K,9,(p) is the forward Alfven conversion coefficient without transition region (see 
Paper I, or equivalently as calculated here with hcor = !)■ 



their exponential asymptotic forms. The two formulae are compared in Figure [Tl indicating that the precise choice 
of coronal structure has little practical significance for the total transmission. Equation (It]) reduces to Equation (44) 
of Leer et al. ( 1982 ) (for reflection coefficient ^ = 1 — T) in the case hcoi- = 1 where a uniform layer is continuously 
appended above an exponential atmosphere. 

Interestingly, even with the hcoi = 100 characteristic of the real solar transition region, substantial Alfven transmit- 
tance (> 10%) is possible if KsecO > 0.1. High frequency, low TR Alfven speed, and large field inclination are all 
favourable for Alfven wave penetration of the transition region. 

However, this is for an incident pure Alfven wave. As shown in Paper I, Alfven waves may also be generated over an 
extended region by conversion from reflecting fast waves, in fact from the evanescent tail beyond the reflection point. 
In this interaction region, they are not pure Alfven waves, only becoming so asymptotically as x — >■ oo. So how do 
these hybrid waves fare on encountering the TR? And what if the fast waves reach the TR before they can reflect, or 
early in their evanescent tail? These issues are addressed next. 

4.2. Alfven Wave Fluxes from Fast Wave Conversion 

The first task is to see how the TR "jump" /icor affects the output fluxes (recall that a^ increases across the TR by the 
factor hcor/h = ft-cor)- Figure [2] shows how the Alfven conversion coeflricient Jl"*" (conversion to the upward propagating 
Alfven wave), reverse coefflcient S^~ (conversion to the downward Alfven wave), reverse fast wave coefficient 0'~ ^ pure 
Alfven transmission coefficient (with TR) 'T and Sl^'T vary with the coronal scale height 1 ^ /icor ^ 100. Here we 
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xr 



= 3.219), XT = 
100, KsecO 



define S^ as the forward Alfven conversion coefRcient without the TR. We fix k = 0.2 (i.e. 
magnetic field incHnation 9 = 30° with wavevector polarizations = 70° and 110°. At h, 
and 'T = 0.033, so a pure Alfven wave is 97% reflected. 

The distance between the reflection point and the TR, measured in units of the chromospheric scale height h 
given by 

a{x^yK^ 



S, and 
0.021 



A 



XT - Xr 



In 



2/l2 



IS 



(8) 



where a(a;^) is the Alfven speed at the base of the TR. This quantity will be seen to be crucial in determining the 
Alfven flux penetrating the TR. To place this dimensionless length in a solar context, let us suppose the chromospheric 
scale height is 150 km, that the density at the base of the TR is 10"^" kgm~^, and that we are dealing with 5 mHz 
waves. Then, A depends on field strength B and dimensionless horizontal wavenumber k = \/(-{(. + 1) h/RQ. It is 
plotted in Figure [3j showing that the fast wave nearly reaches the TR at smaller field strength and spherical degree. 
The latter is due to the wave being more nearly vertical as i gets smaller, meaning that it can reach higher in the 
chromosphere. 
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Figure 3. Distance A = (xt — xji)/h in units of chromospheric scale height between the fast wave reflection height xji and the transition 
region height xt as a function of spherical harmonic degree £ (or dimensionless wavenumber k on the top axis) and magnetic field strength 
B: Full curve 1 kG; dashed curve 100 G; dotted curve 10 G. The wave frequency is 5 mHz, the chromospheric scale height is 150 km, and 
the density at the top of the chromosphere is 10^ ^"^ kgm"''. Where A dips below zero, the fast wave no longer reflects before reaching the 
TR. 



The case /icor = 1 corresponds to the model of Paper I with no transition region. Increasing hcor unsurprisingly 
produces a sharp decrease, to around 1-2%, in the amount of Jl^ propagating through the corona. Interestingly, j^~ 
also decreases, though only modestly. The beneficiary of these decreases is of course !F~ , which increases to over 80% 
at hcor — 100, characteristic of the solar TR, meaning most of the injected fast wave flux returns whence it came as 
a reflected fast wave. Although at (/> = 110° the backward Alfven flux Jl~ is quite significant in the absence of a TR 
(^cor = 1), once /icor rcachcs 100 there is little difference between the two polarizations, consistent with near-symmetry 
of forward and backward Alfven production at = 70° and 110° respectively and the forward Alfven flux having been 
almost totally reflected. 

Also striking in Figure [2] is how well Jl+ is mimicked by ^'T. This appears to suggest that the simple picture of 
Alfven generation according to coefficient ^ followed by TR transmission 'T describes the situation remarkably well. 
This was unanticipated, as the conversion region was not expected to be confined to these few scale heights below the 
TR. 

Next, we fix the magnitude of the transition jump but adjust its distance from the fast reflection point at xr. Figure 
|4] demonstrates the dependence of flux on the location of the TR for flxed k = 0.2, 9 = 30° and /icor = 100 again for 
wavevector polarizations oi (}> — 70° and (/) — 110°. As the TR is raised further above the fast wave reflection point, 
there is a decreasing amount of J^+ with an increase in both downward Alfven and reflected fast wave coefficients. 
The decrease in J^+ is somewhat surprising as we might have expected that a higher xt allows a larger fast-to- Alfven 
conversion region and therefore potentially more outward Alfven flux, at least for — 70° (recall that the conversion 
region is some 20 scale heights wide at k = 0.2). This is apparently because the transition region simply becomes more 
reflective to Alfven waves as the distance between the fast reflection height and the TR increases (see Fig. II]) . Very 
substantial forward Alfven fluxes J1+ of around 20% are seen at small xt ^ Xr = 3.219. It appears that ifthe fast 
wave tail penetrates into the low corona, it continues to generate Alfven waves there. 

One again we see that S^'T closely matches the numerical calculation of J^+. Of course as the transition region 
falls below the fast wave reflection point there are further complications and the two formulae are less well aligned. It 
may be that the Alfven conversion is largely or partially in the corona in these cases, thereby invalidating the simple 
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Figure 4. The forward Alfven conversion coefficient .!?+ (full curve), reverse Alfven conversion coefficient !A~ (dashed curve) and reverse 
fast conversion coefficient ^~ (dotted curve) as a function of xt for n = 0.2, 6 = 30°, /icor = 100, ip = 70° (left panel) and = 110° (right 
panel). The fast wave reflection point is xji = — In 0.2'^ = 3.219 if this is less than xx- Otherwise, the fast wave will propagate throughout 
the chromosphere and immediately be evanescent on entering the corona (provided lj'^ h'^ / a{xj^)'^ < k^). The heavy long-dashed curve 

represents T as given by Eq. JsL The chained lines represent .^1', where J?q (ft, 9, (j>) is the forward Alfven conversion coefficient without 
transition region. 

interpretation of ^'T in terms of two sequential processes. 

We now fix both xt = 8 and hcoi ~ 100, and explore variation with k, 9, and cj). Figure p^ shows graphically how J^+ 
and J^~ vary with transverse wavenumber k and wave polarization (f> for magnetic field inclination 9 — 10°, 40° and 
70°. Across the board, J1+ does not exceed about 2-3%. The backward Alfven flux Jl~ is approximately symmetrical 
about (j) ~ 90° . This is expected in light of the near-symmetry between Jl+ and Sl~ about (j) — 90° found in Paper I 
for the single layer model. We surmise that converted forward-propagating Alfven waves are near-totally reflected off 
the transition region when (j) < 90° and when (j) ^ 90° the fast wave is partially converted to a backward-propagating 
Alfven wave only after fast wave reflection at xu. Maximal backward Alfven flux seems to occur around 50° and 140° 
in each case, but typically at unrealistically high k values {k > 0.5). 

Reducing xt to 2 (Fig. [6| places the TR before (if k < e^"^^/^ = 0.368) or low in the tail of the fast wave reflection 
point, with the effect of substantially increasing A~^ . This suggests that the fast wave evanescent tail, which easily 
penetrates the TR, continues to produce Alfven waves in the low corona. 

Finally, Figure [t] demonstrates how j^+ and R~ vary with n and for fixed 9 — 40° with two different xt values. 
While Jl+ is reduced significantly by increasing xt, there is little change to Jl~. 

4.3. Phase Shifts 

With a transition region in place, we have seen that by far the bulk of the injected fast wave flux is reflected, 
presumably to once again undergo magnetic-to-acoustic mode conversion at the Alfven-acoustic equipartition level 
and then re-enter the Sun's helioseismic wave field. H elioseismic techniques s uch as Time-Distance helioseismology 
( Duvall Jr et al.|[l993 ) and Phase Sensitive Holography ( Braun fc Lindsey|2000 ) use phase to infer travel times and so 
are potentially compromised by any anomalous phase shifts suffered by the fast wave in the chromosphere. The cold 
plasma model affords a particularly convenient testbed for this effect. In such atmospheres the fast wave phase and 
group speeds are just a, independent of the relative orientation of the wavevector and the magnetic field. So, apart 
from the Alfven conversion mechanism, reflected fast wave phase should be independent of 9 and (j). (Specifically, this 
is apparent from Equation (3) of Paper I with ky set to zero. Then the fast and Alfven waves decouple and the fast 
wave is governed simply by (V^ -|- oj'^/a^)£,± = 0, totally independent of 9.) We confirm numerically that indeed the 
phase of the returning fast wave relative to its injection phase is independent of 9 when ky = 0. In this section we 
investigate phase in both the single-layer model and the model with transition region by comparing the phase of the 
returning fast wave with that for the cj) — case where the fast and Alfven waves are entirely decoupled. 

The fast wave is most conveniently characterised by the dilatation, x = ^'^ £^nd thus the forward and backward 
fast wave dilatations x~^ s-nd x~ carry information about their respective phases. The phase change at any point x is 
Sip = — arg( Y"'"Y~). However, we require the difference between the magnetic and non-magnetic cases Aip = Sep — Sipo 
( Cally|[2"009 ), where Sipo refers to the = 0° case, which has no Alfven interaction. This quantity is independent of 
depth provided it is calculated well below the fast wave reflection point eX x = xr. We therefore choose the depth 
where the WKB solution is applied as the point to evaluate the phase difference|_] Figure [s] shows the phase difference 
-180° < Aip s^ 180° as a function of <j) and n for fixed /icor = 100, xt = & and 6* = 10°, 40° and 70°, displaying 



For completeness, we mention the result from Paper I (see also |Ferraro &: Plumpton|l958[ l that the decoupled {ky = 0) fast wave has 
exact solution g^ = a 72^(26"^' ''), with arbitrary complex constant a. This splits into forward and backward travelling wave components: 

- 2 arg a as 



|afl'^^'^'(2e-^/''). Hence x='= = 9x?x = -jafi Ke*''H^^'^'(2e-^'''') + e-'^/^'^H!^^^'!_\{2e-'^/^)]. Consequently, Sipo - 
-oo. This cancels in Aip with an equivalent term in Sip to leave only an effect related to the Alfvenic mode conversion. 
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Figure 5. Left column: The forward Alfven conversion coefficient J?+ as a function of cf> and k for d = 10°, 40°, and 70° (top to bottom) 
as labelled with hcor = 100 throughout. Right column: The reverse Alfven conversion coefficient j?~ for the same cases. In all cases, the 
transition region is at xt = 8. The respective maxima in Si+ are: 0.028 for 6* = 10°; 0.019 for 6* = 40°; 0.023 for 6* = 70°. The bottom two 
panels, displaying the results for d = 70°, are truncated near the top for numerical reasons. This is of little consequence as physically we 
are mainly interested in k < 0.2. 

symmetry about = 90°. The maximum phase difference occurs for higher k as 6 increases. The phase difference is 
generally maximal around = 90°. 

Finally, Fig. [9] demonstrates that for the most part, as expected, phase shift depends only weakly on field inclination 
9. However, this is not the case for cf) « 90°. 

Overall then, it is clear that the phase of the returning fast wave has been very substantially altered by mode 
conversion. Without this interaction with the Alfven wave, A(/3 should be identically zero for all and 0. It most 
certainly is not zero! Any phase sensitive seismology performed using the returning fast wave is therefore severely 
compromised by this highly directional effect. 



5. DISCUSSION AND CONCLUSIONS 
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Figure 6. Same as Fig. Isl but with xy = 2. The horizontal line k = e~^r/2 — 0.368 in the left panels corresponds to the boundary 
between the parameter regions where the fast wave reaches the TR before reflecting (below the line) and where it does not (above). 
Substantially larger values of J?"*", exceeding 0.3, are found below the line, for 6 = 40°. 



We are motivated by doubt about where the Alfven waves that are observed in the chromosphere, corona and solar 
wind are generated. There is reason to believe ( Parkerp991 Vranjes et al. |2008| ) that the creation of these waves in 



the photosphere is problematic. Irrespective of that, it appears from Paper i and Khomenko & Cally (2012) that it is 
easy to generate substantial Alfvenic flux in the chromosphere. The total outgoing Alfven conversion coetticient JV' is 
sensitive to the distance between the fast wave reflection point and the TR. If xt is small, substantial Alfven flux can 
get through, up to about a third, for /icor of 100. Increasing the distance results in stronger reflection of these waves. 

As we see in Figure [8] increasing the magnetic field inclination reduces the impact of the TR on phase shifts. This 
may be due to increasing Alfven path length along the oblique field lines. For the most part, the phase shifts calculated 
here and displayed in Figures [8] and [9] are negative, corresponding to a delay in the phase of the reflected fast wave. 
This could be misinterpreted as an increase in travel time. We interp ret it instead as an art efact of mode conversion. 

How much energy is carried by the observed coronal Alfven waves? Tomczyk et al. ( 2007 ) calculate an insignificant 
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Figure 7. Left column: The forward Alfven conversion coefficient ^+ as a function of and k for fixed = 40° with xj< = 5 and x'y = 10 
(top to bottom) and hear = 100 throughout. Right column: The reverse Alfven conversion coefficient J?~ for the same cases. 



^ though they admit this may be gr oss underestimate due to insufficient spatial resolution. 

(|201l[) using the He ii 304A and Fe ix 



Alfven flux of 0.01 W m 

Recently this figure has been revised dramatically upward by 



Mcintosh et al. 
171 A channels of the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO) to achieve 
arcsecond resolution. They estimate energy fluxes of ^ 100 Wm^^ in active region loops based on typical amplitudes 
of around 20 kms^^. 

In the most ideal case (large attack angle), almost all flux through the a = c level converts to the fast wave. Of 
this, at most 30% is carried through the TR by the Alfven wave, but only if the fast wave reaches or nearly reaches 
the TR before reflecting. Based on the photospheric p-mode power distribution of Thomas ( 1985 ) obtained using the 



Fe I 6303 lin e, which is formed about 290 km above optical depth unity, and the VAL C ernpirical model oflVernazza 



et al 
the 3" 



(|1981|), we crudely estimate {v'^) ~ 4 x 10'^ m^ s~^ and a flux of order 800 Wm~^ associated with p-modes in 
5 niHz band. This appears to be sufficient to supply the coronal Alfven flux estimated by [Mcintosh et al.| 
(|2011|), though the uncertainties in each calculation are large. The link is plausible at least. Of course, one should 



not rule out the possibility of multiple sources of coronal Alfven waves, including direct photospheric excitation and 
chromospheric fast wave conversion. 

APPENDIX 

ON THE UNIDIRECTIONALITY OF PURE ALFVEN WAVES IN AN EXPONENTIAL 

ATMOSPHERE 



Decoupled Alfven waves in the case ky = admit exact solutions i3 



This is because the Alfven wave equation d S,y/dt — a d^y/ds 
into the axisymmetric wave equation on a uniform membrane 



in terms of Hankel functions Hq ' and H^ 



r(i) 



(2) 



(where s is distance along a held line) transforms 






d f d^y 



r dr \ dr 



(Al) 



under the change of variables r — 2ujh/ax, where a^ = a cos 9. The Alfven wave in the exponential atmosphere 
therefore is isomorphic to the problem of axisymmetric 2D waves, with x = +oo mapping to r = 0. As such, there 
is no reflection, despite expectations that the exponential atmosphere might be at least partially reflective to Alfven 
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Figure 8. The phase difference Aip as a function of <p and k for 9 = 10°, 40°, and 70° (top to bottom). Left column: no transition region; 
Right column: transition region is at xj- = 8 with hcor = 100. The 'cut' in the top right panel results from the phases being restricted to 
-180° < Av5 s£ 180°. 



waves. As shown by Whitham ( 1974 S ection 7.4), the 2D axisym metric wave equation with an axial harmonic source 

results in the Hq '{r) solution (see also Courant fc Hilbert|l962 p. 194). We are concerned with the time reverse, or 

equivalently complex conjugate, where x = -|-oo is a sink, selecting Hq (r) instead. See Equation MJ) for the wave 

energy fluxes associated with the two Hankel functions. As a further demonstration that i/g ( r) corresponds entirely 
to wav e propagation in the negative r-direction (positive x), we manipulate Equation (9.1.25) of jAbramowitz fc Stegun 
( 1965 1 to show that 



H, 



(2) 



r)-- 


i 

71- J 


= a{r) + ip{r), 



-r sinh r 



dr 



(A2) 
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clearly demonstrating that there are only Fourier components associated with negative r propagation. The first term 
a(r) = Jo(r) — iHo(r) is a complex Fourier integral containing only waves propagating towards r = 0. The second 
term /?(r) = Ho(r) — Yq^v), is of course real, positive, and monotonic decreasing in r, and so is not wavelike at all. 
In these expressions, Jq and Yq are the Bessel functions of the first and second kind, and Hq is the Struve function. 
Transmission of Alfven waves in an infinite exponential atmosphere is therefore total, T = 1, though discontinuities 
such as we use to model the TR produce substantial reflection (see Equation dsl) 
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Figure 10. The wavelike part a of H^ ' (r) (real part as full curve, imaginary part as dashed curve) and the 'reverberation' part /3 (chained 
curve), all plotted as functions of r = 2exp[— x/2/i]. 



Hollweg & Isenberg ( 2007 1 examined the behaviour of an Alfven pulse in a stratified atmosphere and inferred partial 
reflection in a wake behind the propagating packet. However, we interpret this a little differently. It is well known 
that waves in spaces of even spatial dimension do not obey Huygens' princi ple, meaning that any wav efront trails a 
wake that "persists there indefinitely as a 'reverberation'" in the words of Courant & Hilbert (1962 pp. 208-210). 
Hence, since the Alfven problem is isomorphic to the 2D wave equation, that wake structure is to be expected. 
We see it as represented by /3. It need not be interpreted as a reflection, though of course, one could express /? 

as a Fourier integral, /3 = 2 7r~^ (/rT sin(r cos 6*) dS + J„°° cos(r cosh r) drj representing the reverberation as a sum 

of standi ng waves, consisting of equ al and opposite forward and backward parts. This is consistent with the view 
taken by 
is spatial 



Hollweg & Isenberg (2007). The important point though is that the reflection, if it is interpreted as such, 
y confined to the wake and does not include a finite refiection of energy back to x = — cx). An alternate 



representation of the wake term /3 though is the real integral displayed in Equation (A2), or equivalently the Laplace 
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integral /3 — 2-k^^ J^ e^^'^{l + T^)^^^^ dr. This is an equally valid representation that suggests an in-place oscillation, 
characteristic of solutions of the wave equation in even-dimensional spaces. 

If we turn on a harmonic driver at a particular time, we expect the outermost wavefront to trail a wake which 
develops into the /3 term. Once steady oscillations are set up, after the transient has passed, all energy is propagated 
to infinity. The wavelike and reverberation parts of Hq '{r) are plotted in Figure 10 showing that the latter are most 



significant at small r (large x). Asymptotically, /3 ~ x/hir as a; — >■ -l-oo, so the reverberation takes th e form of a 



straight 'flapping' of magnetic field lines, as seen in the animations accompanying Cally & Hansen (2011). 
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